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Abstract. We analyzed multimode Josephson junctions with exciton-polaritons 
(polaritonic Josephson junctions) when several coupling mechanisms of fundamental 
and excited states are present. The applied method is based on Kcldysh-Green 
'-^J . function formalism and takes into account polariton pseudospin. We found that mean 

Ch I value of circular polarization degree in intrinsic Josephson oscillations and microscopic 

quantum self-trapping follow an oscillator behavior whose renormalizes due to inter- 
mode interactions. The effect of an additional transfer of particles over junction 
barrier occurring in multimode approximation in combination with common Josephson 
■ tunneling is discussed in regime of dynamical separation of two polarizations. 
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1. Introduction 



After theoretical prediction [1] and experimental detection [2] of Josephson tunneling 
between two superconductors separated by a thick insulator under application of an 
external voltage, analogous quantum oscillations were found also in a superfluid system: 
two vessels of superfluid helium connected by a nanoscale aperture [3]. Similar type 
of dynamics was observed for two weakly coupled atomic Bose-Einstein condensates 
(BECs) created in a double-trap potential. Differently from superconductors and 
superfluid helium in Josephson effect for atomic BECs inter-particle interactions play 
essential role [4]. An an-harmonic behavior occurs in these systems additionally to 
common Josephson oscillations. Under certain initial conditions BECs localize in one 
of the traps formed by external potential due to suppression of the tunneling current 
[5]. The phenomenon is known as macroscopic quantum self-trapping (MQST). It is a 
consequence of domination of nonlinearities, e.g. interactions, over Josephson coupling 
in the system. In some other range of the initial conditions delocalized phase is formed 
with unsuppressed oscillations between the traps. 
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In general, such a system - interacting BECs in a Josephson junction (JJ) is 
described by well-known Hamiltonian 

H = H (z(0), 9(0)) = A^f- - cos 0(^1 (1) 

written in terms of population imbalance z(t) = (Nx(t) — iV 1 (t))/(iV 1 (t) + Ni(t)), 
where Ni and N 2 are populations in trap one and trap two, and phase difference 
8(t) between BECs. Two-mode approximation including only lowest energy states of 
symmetric double-trap is used in Hamiltonian (1) [5] and contains single parameter 
A = UqNt/2J. In the last expression J figures as Josephson coupling constant which is 
equal to the difference in energies of symmetric and anti-symmetric states of the double- 
trap potential is less then zero, as the symmetric state is usually lower in energy then 
the anti-symmetric one. Nt is total population of particles in both traps which interact 
with energy Uq. 

Inspecting the phase-space diagram (z,9) of Hamiltonian (1) one can observe two 
distinct regions separated by separatrix line H = 1 ( full/red contour in Fig. 1). Fixing 
the value of H in (1) by choosing initial conditions 2(0) and 0(0), delocalized regime 
will establish if H < 1. It is characterized by the population imbalance being zero in 
average (< z(t) >= 0) for the system's cyclic motion on closed orbits (see Fig. 1). 
Approaching the separatrix from the inner side the orbits starts to deviate from regular 
circles as nonlinearities start to be important. Outside the separatrix where H > 1 the 
population imbalance evolves along open lines with small oscillatins around constant 
mean value < z(t) >= const indicating the transition to MQST effect. 

Besides cold atoms there is an another kind of particles which can undergo a 
transition into collective state with properties similar to Bose- Einstein condensate 
(BEC). These are exciton-polaritons whose condensation has been repported several 
years ago in in CdTe microcavities [6, 7] at about 20 K and later on in GaN-based 
microcavities even at room-temperatures [8, 9]. The possibility to obtain BEC at such 
conditions is extremely interesting from fundamantal point of view. Besides, it opens 
a way for design of novel optoelectronic components, such as polariton lasers [11] and 
various devices based on polariton superfluidity [10] 

Polaritons are two-dimensional quasi-particles appearing due to strong coupling of 
excitons in semiconductor quantum-well(s) and photons confined within a microcavity 
structure. Peculiar properties of exciton-polaritons are consequences of their 
hybrid, half-light half-mater nature coming from photonic and excitonic components 
respectively. In particular, ultra-small effective mass of polaritons, which is typically four 
or five orders of magnitude less then the free electron mass, makes quantum collective 
phenomena very pronounced in polaritonic systems and leads extremely high critical 
temperatures of the polariton BEC transition. 

For the purposes of our further discussion, another property of polaritons is very 
important. That is projection of total angular momentum of a polariton state on 
the structure's growth axis (chosen as z direction) [12]. It can take two values: Jf°' = ±1 
corresponding to mixing of bright excitons (J^ xc = ±1) and photons with right or left 
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Figure 1. Hamiltonian (1) in the phase space (z,9). Contour lines for H=0,l,2 and 
3 arc labeled. Scparatrix is shown with full/red line 

circular polarizations. Dark excitons with J e z xc = ±2 are uncoupled from cavity mode 
due to optical selection rules and thus do not contribute to the formation of a polariton 
doublet. From the point of view os spin polariton is thus analogical to electron, as both 
of them are two-level systems. Consequently one can apply pseudospin formalism for 
description of spin dynamics of the polaritons. It is convenient to represent a pseudo 
spin state of polariton system by a point on Poincare sphere (also known as Bloch 
sphere). Importantly, polariton pseudospin unambiguously defines the polarization state 
of photoemission from the cavity. The states lying at the poles of Poincare sphere will 
correspond to circular polarization of the polaritons, the states on equator- to linear 
polarizations and all other states to elliptical polarizations. 

Being bosons, exciton-polaritons should manifest phenomena related to Josephson 
effect. Polaritonic Josephson junctions (PJJ) were considered theoretically in Refs. 
[17, 18, 19] and experimentally in Ref.[22]. In these articles polariton pseudospin 
was neglected, and and Josephson coupling occurred between the polaritons located 
in spatially separated traps. The introduction of polarization degree of freedom makes 
Josephson dynamics far more rich. Indeed, in addition to the coupling between states 
localized in different traps, two polarization states in single asymmetric trap can be 
coupled as well by mechanism analogical to TE- TM splitting [21]. Josephson-type 
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oscillations can occur between different polarizations, the phenomenon which was called 
intristic Josephson effect [20]. The interplay between intrinsic and extrinsic Josephson 
dynamics can lead to dynamic separation of different polarizations in the real space and 
other intriguing phenomena [20]. 

The role of the polarization coupling in PJJ was studied for several configurations 
[23, 24, 25]. In most of them so called two-mode approximation was applied, in which 
only two lowest levels, one in each of two traps, were considered. An exception is 
the reference [25] where several modes were included in analysis in order damping of 
Josephson oscillations due to the interactions of the polaritons with acoustic phonons. 
However, even in this case the separation between two fundamental modes and excited 
ones was considered to be much greater then characteristic values of blueshifts provided 
by polariton- polariton interactions, and transitions to the excited modes were only 
possible due to absorbtion of phonons. 

In this paper we analyze multimode PJJ with confining potential created in such 
a way that interactions between the lowest and excited levels are not negligible. This 
kind of coupling have been already considered for Josephson Junctions based on atomic 
BECs [13]. Our goal is to extend the analysis of Ref. [13] for the case of cavity polaritons 
accounting the polarization degree of freedom, and clarify the effects of the coupling 
between fundamental and excited modes on nonlinear polarization dynamics of PJJs. 

2. Model 

In Ref. [13] Keldysh-Green functions technique was employed to study dynamics for 
atomic Josephson Junction with multimode structure. Here we shall follow similar 
approach for description of PJJs accounting for the spin degree of freedom. Because of 
spinor nature of polariton condensate in addition to normal Josephson coupling we have 
to consider also and the coupling of two polarization components, normally irrelevant 
for atomic JJs. 

We start by expanding exciton-polariton field operators ^f a ± (r, t) over complete set 
of eigenstates {fx, f2, ■■■} of a double- well potential 

oo 

*a±(r,t) = h(r)a la± (t) + f 2 (r)a 2a ±(t) + ]T f n {r)b nfT ±{t) (2) 

n=3 

where a la ± and a 2cr ± denote annihilation operators for two lowest modes of PJJ localized 
in the traps one and two, and b rm ± stands for the anihilation operator on excited level 
n taking values, denotes circular polarization of the state. 

The system we study is described by the Hamiltonian which can be represented as 
a sum of three terms: 

H = H a + + H a - + H a + a - . (3) 

The terms conserving z- projection of pseudospin a + and <r_ denoted by H a + and 
H a - are given by following expressions 

H a± = J d\¥ a± {r,t){-^ + V ext {r,t))i> a± {r,t) (4) 
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+-j d 2 r J d 2 r'q+ ± (r,t)¥ ± (r\t)V{r - r')£ CT ±(r' ,t)$ a ±(r,t). 

The effective mass of polaroton is denoted by m and V ex t(r, t) stands for external double- 
well potential. V(r,r') = g5(r — r ) is contact interaction described by delta function 
and interaction constant can be estimated as g ~ Ebcl 2 b with Eq and a B being being 
exciton binding energy and Bohr radius of the exciton respectively [14]. 

The last term in the Eq. (3) accounts for the Josephson- type coupling of the 
particles having opposite pseudospin projections. It can be viewed as a consequence of 
one presence of effective in-plane magnetic field Q(r) arising from the asymmetry of the 
structure [21] and acting on polaritons spins and can be represented as 

iW = / d 2 r¥ ± (r,t)Q(r)q^(r,t). (5) 

After subtitution of the expensions (2) into the Hamiltonians (4) and (5) the starting, 
i.e. full Hamiltonian (3) can be recast as 

H = H° + H exc + H int . (6) 

where the first term describes the dynamics of the four fundamental modes (accounting 
to polarization degreee of freedom), the second term describes the dynamics of the 
delocalized excited modes and the last term corresponds to the coupling between 
fundamantal and excited modes. These terms read: 

H° = E Y^ a\ a a i<T + J ^(a^^ + a^a^) + Y a\ a a\ a a ia a i(J + fi Y aL^i-a! (7) 
H exc = £ (E n + U nm < blJ ma >)b{Jb rm + Q n Y, biX-* (8) 

n,m;<j n;cr 



run 
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[< b lnab}na > K0K0+ < b ma b ma > fe^L); 



H mt = K n J2 [\ (a\AXX. + h.c.) + 2a\ IJ a i XXA + (9) 



i,n;a 



nob no 

+ h.C.)}. 

n;o 

In the above expressions E n are energies of the modes, Jo is a Josephson 
coupling strength between two fundamental modes localized in right and left traps 
of the double- well potential, fl n are coupling strengthes between states of different 
polarizations at level n corresponding to intrinsic Josephson effect, Uq describes 
polariton- polariton interaction in the fundamental modes, U mn - interactions in excited 
modes, K n - interactions between fundamental and excited modes and J n describes the 
renormalization of the Josephson tunneling due to the interaction between fundamental 
and excited modes. The parameters entering into Hamiltonian 6 are related to those 
entering into the initial Hamiltonian 4 as 

E oi n) = J d 3 rfl 2{n) (r) f-|-A + V ext (r)) /i, 2( „)(r), (10) 



Jo = J d 3 rfl 2 (r) (^-A-A + V ext (r)j f 2jl (r), 
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(11) 

U = J d 3 rdVfl 2 (r)fl 2 (r')V(r,r')f 1 , 2 (r')f 1 , 2) (r) = 

= 9 I d 3 rf m (r)\ (12) 



U rm = gJ d 3 rf£(r)f*(r) } (13) 
^o(n) =9 J d 3 rfl^ 2a(na) (r)n{r)f\_^ 2 _ a{:n _^{r) } (14) 
K n = 9 f d 3 rf* 2 2 (r)f n (r) = g [ <Pr\fo(r)f n (r)\ 2 , (15) 



J n = 9 J d 3 rf*(r)f 2 (r) |/ n (r)| 2 = g J d 3 rf*(r)f;(r)f n (r). (16) 

In the part corresponding to the excited states H exc the interactions are treated 
using the mean field approximation [26] 

In order to obtain a closed system of dynamical equations describing interacting 
PJJ we write Keldysh propagators of the system in the following representation 

(t,t) G a l3- a -cr(t,t ) J 

where the general indices a and /3 become i or j for the ground states and take values 1 
or 2 as there are two lowest lying modes. The excited states are counted by associating m 
or n to both a and (3. As previously indices ±er denotes pseudospin degrees of freedom. 
The elements of the above matrix (18) are themselves 2x2 block-matrices of the form 

■A ( + /\ | G a p aa (t,t) F a pcrcr(t,t) \ , v 

iGapaaKht ) = \ p , ~ , r . (19) 

with the elements being, for example for the reservoir modes 

GWM) =< ^b la (t)blAt) > (20) 
F a p aa (t,t) =< ^b l(T (t)b ma (t) > . (21) 



and similarly for the fundamental modes localized on the traps 1 and 2. Time-ordering 
ST in the previous formulae is performed on the Keldysh contour [27] and appears 
because of non-adiabatic switching of Josephson coupling in the initial time instant. 
As it is not possible to guarantee the behavior of the system with such a kind of the 
irreversibility for asymptotically large times [28] the time contour in Keldysh formulation 
is adapted so that system evolves forwardly on time axis to some time and then from 
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this point it makes backward evolution to the initial state. Any Keldysh time-ordered 
product of two operators, for example, operators b n and b m , has form 

< &[b m (t)b n (t')] >= 6(t,t') < b m (t)b n (t) > +6(t',t) < k(t)b m (t) > -(22) 

where 9(t,t') is "step" function defined on the Keldysh contour for two arbitrary times 
t and t in way that is one always when the first argument is later then the second one, 
and zero otherwise. In this sense the first Green function in the formula (22) is called 
"greater", denoted usually with F > and the second one is "lesser" Green function F < . 
They act only on the forward or the backward Keldysh contour branch, respectively. 

Using Wigner transformation, the matrix elements G a p aa (t,t') can be written in 
terms of center of "mass" and "relative" time coordinates T = (t + 1')/2 and r = t — t' . 
We are interested here in external and internal Josephson dynamics which are much 
slower processes then the others occuring in the system so that we can work in the limit 
in which r = 0. The new matrix elements G a p aa {T,r) will then only depend on the 
macroscopical time T [13]. According to the expression (22) we will then deal only with 
"lesser" functions F<(T) and G<(T). 

Equations of motion techniques combined with use of the mean- field approximation 
allows us to obtain the closed system of equations for Green functions defined by Eq.19: 



r\n< 

' ^rp ^na r nana ^na r nana ! > c "n\^nan-a ^n—anaji \^°) 

i^~Pr ^na ~ '^n-a)G nan _ a = 2^ na F nan _ a — 2^ na F Tmn _ a — f2 n (G — G n _ an _ a )] (24) 



^rp 1 -«» -ii—vj — nan— a nan— a nan— a --'<-\—najc 

d_ 

dT 
d 

n—an—a > 



jrn 2E na 2T no .) F nana — 2^ na G nana + 2^ na G nan(7 Vl n F^ an _ a ; (25) 



(i ^j-, + 2E na + T ncr + r T n ^ rj )F nan _ a — 2[p, na c n — a )G rum _ l7 Q, n {F n(Jna + F 1 



(26) 



where n — 3,4, ... stands for the n th excited level. Dynamics of the fundamental states 
i = 1,2 is given by following expressions: 

i* — (-^o + U Gf aia . + 2A n 2J Gnana)Gf aiu + ( J + 2 J n ^ Gn anr7 )Gf aja + (27) 

na na 

+i(K n F iaia + J n F- a - a ) ^ F nana ~ QoG iai _ a - 

na 

with i,j = l,i ^ j. 

Self-energies T and H read 

-zT n(7 = U*J2 G< ama + 2K n (G< + G 2 <) + 2 J n (G< 2a + G< J, (28) 

m 

— i^na — — X] ^mama ^ ^"(^l^-lo- + ^2a2a) + ^n^Xola' (29) 



The normal and anomalous self-energies represent the energy renormalizations entering 
to the diagonal and off-diagonal propagators in the matrix (19) due to particle- particle 
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interactions present in PJJ. For simplicity we take them diagonal in excited level 
index n neglecting collisions between particles situated at different reservoir levels, 

v nm v u nm . 

The system of (23)-(29) is the closed system of nonlinear first order ordinary 
differential equations which can be solved numerically. The corresponding analysis is 
presented in the next section. 



3. Results and Discussions 



We consider a PJJ inside GaAs- based quantum microcavity with parameters similar to 
those used in Ref.[20]. We numerically studied the system of equations (23)- (29) in order 
to analyze the effects of the interactions between fundamental and excited states given 
by parameters K n and J n on various types of Josephson deynamics. In the reference [13] 
it was found that for atomic condensates such kinds of interactions in general lead to 
chaotization of the Josephson oscillations after some initial period of regular dynamics. 
The effect is due to the intensive exchange of the particles between fundamental states 
and multi- mode reservoir of the excited states. Here we consider time intervals smaller 
then those necessary for the transition to chaotic regime. The reason is that polaritons 
have finite lifetimes, and in the regime of pulsed excitation they will simply disappear 
before the system will demonstrate the characteristics of chaos [20]. 
The quantity 

(rp\ ^*na+ia+ CO ^n— ct~CO ( c i(\\ 

Pn{-L ) - T7< nf\ JT< 77fW l JU J 

^na+ia+l 1 ) ~ ^n-a-i-o-K 1 ) 

describes circular polarization degree at the state n. Its dynamics is shown at Figure 
2. For a while, we neglected the extrinsic Josephson coupling puting J = 0, J n = 0. 
Panels a) and b) show a profile of the oscillations of circular polarization degree in 
the fundamental states corresponding to intrinsic Josephson effect and corresponding 
Fourier spectrum. The dashed line correspond to the case when coupling to excited 
states is switched off, while solid line accounts for this coupling. One sees, that 
introduction of the term K n slightly renormalizes the frequency of the oscillations. 
However, the effect remains quite weak, as population of the excited levels is sufficiently 
small, less then ten percents of total number of particles. On the contrary, the effect 
of coupling becomes more pronounced if one monitors the intrinsic Josephson dynamics 
on the excited states, shown at panels c) and d). The interaction with fundamental 
modes changes the oscillation pattern on excited modes quite radically. The effect is 
clearly seen at panel d) showing Fourier power spectrum of the oscillations. Account 
for the terms K n leads to the appearance of higher harmonics in the spectrum. Besides, 
in place of the intrinsic Josephson oscillations with zero time average, novel regime 
establishes in which < p(T) >^ 0. In the absence of coupling K n the frequencies of the 
intrinsic oscillations in the ground state and the reservoir are different (dashed lines in 
panels a) and c)) although polarization couplings on these levels are equal fl = Q n . It 
is result of the renormalization of the oscillation frequency for the fundamental mode 
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Figure 2. Panels a) and b) Profile of the intrinsic Josephson oscillations of circular 
polarization degree at fundamental states and its Fourier power spectrum. Panels c) 
and d) Profile of the intrinsic Josephson oscillations of circular polarization degree at 
first excited state and its Fourier power spectrum. Dashed lines correspond to the 
absence of the coupling between fundamantal modes and excited states, K n = 0. Solid 
lines correspond to the case K n = 0.3Uq- 

due to nonlinearities, which is negligible in the low populated reservoir. Account of 
the reservoir- fundamental mode coupling makes the frequencies of intrinsic Josephson 
oscillations comparable for all modes of the system. 

Fig. 3 illustrates behavior of the polarization degree when ground state MQST 
phase in the intrinsic Josephson effect interacts with the reservoir. As compare the 
situation analyzed in Fig 2 particle-particle interactions characterized by a parameter 
Uq dominate over Josephson coupling Q in the fundamental mode . The profiles of 
the oscillations of circular polarization degree for the fundamental modes are shown in 
the top of the panel a). The solid line corresponds to the case K n ^ and dashed 
one for K n = 0. The dashed line is characterized by a single peak in the Fourier 
spectrum, usual for MQST regime in two- mode PJJ. On the contrary, the introduction 
of the coupling with excited level leads to the appearance of the two additional peaks 
in Fourier spectrum as it is shown at the panel b). The similar trends can be seen for 
oscillations in the excited level, illustrated by lower curves at panel a) and blue line at 
panel b). 
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Figure 3. Top of the panel a) shows the polarization degree of the intrinsic Josephson 
oscillations in the fundamental mode. Bottom of the panel a) shows the same quantity 
in the reservoir. Dashed lines: K n = 0. Thick/blue and solid/red lines stands 
for K3 = O.SUq. Panel b) shows corresponding Fourier power spectrum when the 
interactions are present (red for the fundamental mode, blue for the reservoir). 



Finally we consider the scenario in which two fundamental modes localized in traps 
one and two, initially elliptically polarized, are coupled with each other by means of 
extrinsic Josephson tunneling J. It was already pointed out that spatial separation 
of polarization occurs in this case [20] for sufficiently high polariton concentrations 
and circular polarization degrees. In this regime the dominant polarization component 
is trapped in one of the wells (dotted line in lower panel of Fig. 4), while another 
one undergoes Josephson oscillations (dotted line in upper panel of Fig. 4). In the 
presence of the higher modes there is an extra pseudospin-conserving exchange of the 
particles through the excited levels between the traps described by reservoir- assisted 
coupling term J n in Eq.9). The influence of this term is illustrated in Fig. 4 by solid 
red and blue lines corresponding to J n = O.IUq and J n = 0.3Uo respectively. For freely 
oscillating component the presence of the reservoir-assisted term leads to self-trapping. 
This transition is provided by competition of two effects: common Josephson tunneling 
J and reservoir- assisted coupling given by J n . As their signs are opposite (J < and 
J n > 0), increase of J n decrease in fact absolute value of the total Josephson coupling. 
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Figure 4. Upper panel: Josephson oscillations of the population imbalances between 
two traps for minor polarization component. Lower panel: Josephson oscillations 
of the population imbalances between two traps for major polarization component. 
Jo = 50fieV. Dashed line: no reservoir assisted coupling, J n = 0. Red line: J3 = O.IUq, 
blue line J 3 = 0.2U . 

This results in overwhelming of free oscillations by nonlinearities and establishing of 
MQST in the minor polarization component. The smaller is the total Josephson coupling 
the trapping becomes stronger (compare blue and red line of the upper panel). The 
dominant component rests relatively robust to the influence of the reservoir-assisted 
coupling. Note, that in hypothetical case of positive J the result will be opposite to 
those considered here: the reservoir- assisting terms will increase the absolute value of 
the effective tunneling constant, thus contributing to the distruction of MQST regime 
for major polarization component. 

4. Conclusions 

In conclusion, we analyzed polarization dynamics in multimode Polartonic Josephson 
Junctions. We have found that the coupling between fundamental and excited modes 
changes the patterns of intrinsic and extrinsic Josephson oscillations, leading to the 
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appearance of higher harmonics and transitions between Josephson and MQST regimes. 
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